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^ i Information about the physical processes that generate the primordial 

'-^1 fluctuations in the early universe can be gained by testing the Gaussian 

I ' nature of the fluctuations through cosmic microwave background radiation 

2 ■ (CBR) temperature anisotropy experiments. One of the crucial aspects of 

^ . density perturbations that are produced by the standard inflation scenario 

is that they are Gaussian, whereas seeds produced by topological defects left 
over from an early cosmic phase transition tend to be non-Gaussian. To 
^ ' carry out this test, sophisticated statistical tools are required. In this paper, 

■ we will discuss several such statistical tools, including multivariant skewness 

and kurtosis, Euler-Poincare characteristics, the three point temperature cor- 
relation function, and the Hotelling's statistic defined through bispectral 
estimates of a one dimensional dataset. The effect of noise present in the cur- 
rent data is discussed in detail and the COBE 53 GHz dataset is analyzed. 
Our analysis shows that, on the large angular scale to which COBE is sensi- 
tive, the statistics are probably Gaussian. On the small angular scales, the 
importance of Hotelling's statistic is stressed, and the minimum sample 
size required to test Gaussianity is estimated. Although the current dataset 
available from various experiments at half-degree scales is still too small, im- 
provement of the dataset by roughly a factor of two will be enough to test 
the Gaussianity statistically. On the arcminute scales, we analyze the recent 
RING data through bispectral analysis, and the result indicates possible de- 
viation from Gaussianity. Effects of point sources are also discussed. It is 
pointed out that the Gaussianity problem can be resolved in near future by 
ground-based or balloon-borne experiments. 



^Address after Oct. 1, 1993 will be: 301 LeConte Hall; Center for Particle Astrophysics; 
University of California; Berkeley, CA 94720. 



1 Introduction 



The cosmic structure formation problem is essentially an initial value prob- 
lem: how did the universe generate the initial perturbations? In particular, 
one can divide the initial condition models into two clear classes: Gaussian 
or non-Gaussian? Cosmic inflation |l|], on one hand, provides a natural way 
to generate Gaussian initial perturbations |0] ; spontaneous symmetry break- 
ing, on the other hand, will lead to the formation of topological defects 
via Kibble mechanism and the perturbations generated by topological 
defects can be characterized as non-Gaussian. Non-Gaussian perturbations 
also arise in various non-standard inflationary models ||^. Thus, a test of the 
Gaussian nature of the primordial perturbations will not only be helpful in 
discriminating different models for structure formation, but could also shed 
light on new physics that yield topological defects or special non-standard 
inflations in the early universe. Such a test is, therefore, very important and 
timely in today's cosmology. 

There are two ways to carry out the test. One is from the statistics of 
the galaxy counts in a redshift survey ^ . However, since the density fleld we 
observe today has already gone through the "black-box" of non-linear grav- 
itational evolution, one has to fllter out this effect carefully to get a reliable 
estimate of primordial quantities [0]. In this paper, we will concentrate on 
the other approach, which is from the cosmic microwave background radia- 
tion (GBR) anisotropy experiments. The experiments measure the primeval 
density perturbations at redshift z ~ 1000. The density contrast is fairly 
small at this epoch. The Gaussian nature of the microwave background fluc- 
tuation directly reflects the nature of the primordial perturbations. This 
approach is promising, especially after COBE's detection of the temper- 
ature anisotropics at large angular scales, and the continuing accumulation 
of data on smaller angular scales 

Prior to COBE's detection, studies on CBR were focused on determining 
the level of anisotropics. The Gaussianity test of the anisotropics was largely 
considered as a next-step problem and experimentally intractable. Few de- 
tailed studies [jlO| [|ll|] on the statistics of the CBR anisotropics have been 



carried out except for the Gaussian case |12|. Now that the CBR anisotropics 



are detected [||||T3|) this next-step but important question should be brought 
into focus and we are optimistic that it can be resolved experimentally in 
near future. As we will show later, one doesn't need a full-sky coverage at 
high angular-resolution to study the Gaussianity problem statistically. 

There are several recent papers on statistical tests for Gaussianity fl^ 
To carry out the test through CBR experiments, two important feature have 
to be stressed. One is the instrumental noise which is present in all current 
experiments. One should have a clear understanding of the instruments and 
associated noise before attempting to decipher non-Gaussian signals from 
experiments. The other feature is the smoothing scale 6s each experiment 



operates. Note that the angular size of the comoving horizon at decouphng 



epoch is 9c ~ 2° by assuming standard recombination |Tj]. Put this charac- 
teristic scale in mind, one can divide all CBR experiments into three cate- 
gories: large scales {6^ ^ 6'c), intermediate scales {Og ~ 1°) and small scales 
{9s ~ arcminutes). For large scale experiments each measurement is 

a sum of anisotropies in several independent horizons, and one would expect 
the statistic to be close to a Gaussian simply by the virtue of central limit 
theorem . For small scale experiments , we will show later that the 

data fails the Gaussian statistical tests. However, on arcminute scales, fore- 
ground source contaminations are important. The statistics of the data may 
not reflect the statistics of the CBR anisotropies at these scales. Interme- 
diate scale experiments pi| are ideal for testing Gaussianity. Although the 
dataset available is still too small, as we will show in section 4, improvement 
of the dataset by roughly a factor of two will be able to test the Gaussinity 
of CBR statistically. 

Several statistical tools are discussed in this paper. In section 2, we dis- 
cuss the simplest tests of Gaussianity through the skewness and kurtosis of 
the one point distribution. Skewness and kurtosis are the normalized third 
and fourth moments of the distribution and they vanish for Gaussian distri- 
bution. Several physically motivated non-Gaussian probability distribution 
functions (PDFs) are considered, and effects of noise are discussed. To con- 
sider the possible correlation between the signal and noise, multi-variant 
skewness and kurtosis are introduced and their statistics are discussed. In 
section 3, one geometrical measure of the random field, the Euler-Poincare 
(EP) characteristic, is discussed and used to test Gaussianity. The statis- 
tics of the EP-characteristic and the effects of noise are discussed and we 
show that EP-characteristic is hardly a good discriminator between Gaus- 
sian and non-Gaussian fields when the noise is comparable to the signal. 
In section 4, we discuss using the three point correlation to test Gaussian- 
ity. Theoretical predictions in various models are discussed and we present 
our analysis of COBE 53 Ghz data. The result is in good agreement with 
Gaussian assumption. In section 5, we discuss using Hotelling's statistic 
to test Gaussianity on intermediate angular scales {9 ~ 1°). The minimum 
sample size to carry out the test is estimated and sampling technique are 
also discussed. Although the current dataset is still too small to carry out 
the Gaussianity test, improvement of the sample size by roughly a factor of 
two will be adequate. In section 6, we use the statistic to test RING 
dataset from OVRO ||2^ on small angular scales (~ arcminutes). It is found 
that the data is not consistent with Gaussian distribution. However, one 
cannot conclude that the CBR anisotropies are non-Gaussian on these scales 
because of the foreground source contaminations. In section 7, we discuss 
looking for a special non-Gaussian signal, the point-like CBR anisotropy, in 
small scale CBR experiments. The Gaussian nature of perturbations from 
inflation is shown in the appendix. 



2 Skewness and Kurtosis of Noisy Data 



The simplest tests of Gaussianity will be skewness /is and kurtosis /X4 of 
the distribution of temperature anisotropies S, 

/i3 = "^3/o-^, /i4 = m4/cr^ - 3, (1) 

where and 7714 are the third and fourth moments of the distribution, adn a 
is the variance of S. For Gaussian distribution, both /is and fi^ vanish. In this 
section, we discuss several physically motivated non-Gaussian distributions: 
exponential, log-normal and x^- As we expect, noise will blur the effects of 
non-Gaussian distribution. The skewness and kurtosis for these distributions 
are calculated both with and without noise. We also discuss the use of 
multivariant skewness and kurtosis in cases of noisy data and show how to 
estimate them from experimental data. 



2.1 Skewness and Kurtosis of Non-Gaussian Signals 



Coles and Barrow ||T0[ have studied the statistics of a large class of non- 
Gaussian distributions. We choose the following distributions based on phys- 
ical considerations. To reflect the real experimental setup where the mean 
of the distribution is subtracted, we standardize the distribution so that all 
of them have a zero mean. Furthermore, we normalize the variance of the 
distribution to be unity. Thus, distributions we study below correspond to 
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the probability distribution functions (PDF) oi x = ^ ■ —, where cxo is the 



observed rms temperature fluctuations. 

1. Exponential distribution. 

This distribution may describe the temperature fluctuation produced 



from the cosmic string network on arcminute scales |Tl||. The PDF of this 
distribution is; 

P(x) = -^exp(-v^|a;|). (2) 
The skewness and kurtosis of the distribution are: 



/is = 0; /i4 = 1.5. (3) 

2. Log-normal distribution. 

This distribution is widely used in the statistical studies of galaxies and 
clusters. It is temping to suggest that it might describe the distribution 
of temperature anisotropy from Sunyaev-Zeldovich (SZ) effects Since 



the effect is produced by the hot-gas in the rich clusters, thus it should 
relate intrinsically to the distribution of rich clusters, which is log-normal. 
Simulation of SZ in the cold dark matter scenario seems to support this 



connection. A small reminder is that the SZ effect always produces cold spots 



on the sky; thus, the distribution of the temperature fluctuation is different 
from the usual log-normal distribution by a sign. The PDF is given by: 



^(^) = , . exp(-(Zo(7(-x))V(T2),x< 0. (4) 



where a is given by: 

. = ll„i±^, (6) 

2 2' ^ ^ 

The skewness and kurtosis of the distribution are: 

exp(3a) — 3exp(a) + 2 



(exp(a) - 1)3/2 



exp(6Q;) — 4exp(3a) + 6exp(Q;) — 3 
= (exp(a) - ly ' 

where a = cr^. Thus, 

/i3 = -0.66, /i4 = 2.72. (7) 

3. Xn distribution. 

This class of distributions provides a good fit to the statistics of tempera- 
ture fluctuation from global topological defects and non-topological defects in 
the framework of 0{N) cr-model [^. In this model, a global symmetry 0{N) 
is broken to 0{N — 1) by a N-component real scalar field = (01, 0Ar) in 
the early universe. The temperature anisotropy produced by the dynamics 
of the scalar fields is given by: 

Y = — = — ^ Y.[<P^ + (8) 

When is small, the dynamics of the scalar fields are nonlinear, thus, the 
PDF for (pi is non-Gaussian. However, when is larger, the dynamical equa- 
tions for (pi decoupled and become linear ||2^ and therefore (pi{i = l,...,n) 
become independent Gaussians. From Eq. (H), it is clear that in the large N 
limit, the temperature anisotropy is distributed with n = 4N degrees of 
freedom. Because of our standardization process, the PDF is related to the 
usual distribution by the following transformation: 



+ a'^^l/V2n. (9) 

Thus, the PDF of the distribution is given by: 



/X 1 / / — T\ 2-1 / V 2n ■ X + n , / — ~ 
= (2/n)W2r(n/2) ^^+V"/'^ ' ^"^^ 2 ^'^ " 

(10) 

This PDF is plotted in Fig. (1) for n = 4, 8, 16. The distribution is very 
non-Gaussian for low n but approaching asymptotically toward Gaussian in 



the large n limit. The skewness and kurtosis of the distribution of this class 
of distribution is: 

/is = -y/s/n, /i4 = 15/n. (11) 

With the presence of Gaussian noise, the skewness and kurtosis will reduce 
dramatically. While the third and fourth moments of the distribution are 
unchanged, the variance increases by (1 + q;^), when the signal to noise ratio 
is 1 : a. Thus, the skewness will reduce by (1 + a^)"^/^ and the kurtosis 
by (1 + For the log-normal distribution, the relations mentioned above 
are not true. The skewness and kurtosis can be calculated from Eq.(^. The 
skewness and kurtosis for various non-Gaussian distributions are listed in 
Table 1. In the noisy case, the signal to noise ratio is set to be 1:1. 

2.2 Estimate of Skewness and Kurtosis of Noisy Data 

In this section, we will address the following questions: how to estimate 
the skewness and kurtosis from noisy experimental data, and what are the 
statistics of these quantities when the signal and noise are both Gaussian? 

In the case where the experimental data is noisy, one has to deal with 
two random variables, the signal and the noise, which have a bivariant joint 
distribution. Thus, we have to generalize the usual skewness and kurto- 
sis to bivariant distributions. Let us first consider the general multivariant 
distribution. For a p— dimensional random vector X = {xi,X2, ■■■,Xp) with 
zero mean and covariance matrix S, it is helpful to introduce the following 
multivariant measures of skewness and kurtosis [^: For skewness, 

p p , , , 

n rr ^ss tt ,,rst , ,r s t /-in\ 

P^p= 1^ 1^ ^ ^ ^ /^iii/^iii > (12) 

where /i"* =< XrXsXt >, cr*-' is the i,jth element of the inverse of 

covariant matrix S. If the joint-distribution function for X is a multivariant 
Gaussian, then (3ip = 0. (3ip can be estimated from a sample of size N, where 
we can replace the ensemble- averaged /im and covariance matrix elements 
with the sample-averaged onesQ. If we denote the estimate of Pip by bip, 
then under the Gaussian hypothesis, the statistic A = is approximated 
distributed as a x 

with EiP±MP±Il de grees of freedom. For us, the interesting 
cases are for p = 1,2. When p = 1, = fj,^, the usual skewness. If one 
denote 63 as the estimate of skewness fi^ from a sample of size A^, then 

b^ = ^hl (13) 

i=l 

^Here, we explicitly assume that the temperature fluctuation is ergodic so that the 
ensemble average will be identical to the spatial average (in our case, it is the sample 
average) in the limit of large spatial coverage. The ergodicity is guaranteed if the power 
spectrum of the fluctuation is continuous. See Adler and Bardeen, Bond, Kaiser & 
Szalay ^ for details. 



and A = is distributed as a with 1 degree of freedom. For p = 2, 

(3u = (1 - pTHiIo + 7o3 + 3(1 + 2p'){^l, + 72^) - 2p3^3o7o3 

+6p[(73oP7i2 - 72i) + 703(P72i - 712) - (2 + p^)7i272i]}, (14) 

where 

P = ,lrs = ^T^, Prs =< X^X^ > ■ (15) 

0-10-2 0-[0-| 

In the case where signal and noise are both random Gaussian fields, the joint 
probability distribution function is a bi-variant Gaussian: 

^^^''^'^ " {2iTYdetM ^xp(-V2E^»^ii'^i)> (16) 

where the correlation matrix between the signal and noise is the following: 

Mn = al, M22 = al, = M^i =< X1X2 > . (17) 

For COBE DMR, the correlation matrix can be estimated from (A+B)/2 and 
(A-B)/2 maps. The bi-variant measure of skewness jSu = if we assume that 
the signal and noise are both Gaussian. Let us denote 612 as the estimate 
of (3i2 from a sample of size (say, COBE map). Under the Gaussian 
hypothesis, A = ^^^^ is distributed as a with 4 degrees of freedom. 
For kurtosis, we have the following measure: 

/?2p = [< XS-^X >]2 - p{p + 2). (18) 

Given a random sample of size N on the random vector X, we can replace 
the ensemble average with the sample average and estimate (322 by 

I N _ 

K = ^ - p)S~\X, - ^^)f - p{p + 2), (19) 

where p = jj Yli Xj, and Sij = Yli {Xi — p){Xj — fj). For a multivariant 
Gaussian, h2p = 0, and the following statistic B = [ sj(pl2) V^"^ is approximately 

JV 

distributed as normal with zero mean and variance one. 

For p=l, P21 reduce to /i4, the usual kurtosis we defined before. An 
estimate of fi^ from a given sample is 

1 ^ 

&4 = ^E^'- (20) 

i=l 



and the sample variance is \J2A/N . 
For p =2, 

(322 = -8 + {740 + 7o4 + 4p(p722 - 713731) }/(l - P^)^ (21) 

where p and 7^3 is the same as the ones defined in Eq. ( |T5|). Under Gaussian 
assumption, (322 = 0, and the estimation of /?22, ^22 has zero mean and sample 
variance ^JaS/N. 



2.3 COBE data 



At present, the only existing complete dataset of temperature anisotropies 
is the COBE dataset. Thus, it is tempting to ultilize the statistics we dis- 
cussed before and put a constraint on the possible deviations from Gaus- 
sianity through COBE data. Unfortunately, the skewness and kurtosis of 
the distribution of ^ cannot be estimated directly from the COBE dataset 
because the signal in each pixel is not an independent measurement of 
Furthermore, in order to estimate skewness and kurtosis, we have to assume 
that the temperature anisotropies are ergodic so that the spatial average is 
equivalent to the ensemble average. But, even for Gaussian fluctuations, the 
ergodic hypothesis is true only if the power spectrum of the fluctuations is 
continuous The observed temperature anisotropies are 2D random 

fields on a 2-sphere. The power spectra Ci is discretized and asymptotically 
approaching continuity in the large I limit. For the COBE experiment, which 
is sensitive only to the low / moments, the ergodic hypothesis doesn't hold. 
Thus, the statistics estimated from COBE directly will only be the estimates 
of the "local" values: they are the measure of deviation from Gaussiantiy in 
our horizon. To estimate the cosmic skewness and kurtosis, i.e., the skew- 
ness and kurtosis averaged over an ensemble of horizons, one has to treat 
the observed value as a N-dimensional random vector, where N is the sample 
size (the total number of pixels). Since we have only one measurement (we 
have only one universe), in order to calculate the ensemble average of the 
quantities equation, the only conceivable way is to use Monte-Carlo simula- 
tion. The COBE analysis along this line will soon appear |^ and won't be 
repeated here. We note that by using the statistics of multivariant skewness, 
which is a distribution with A^(A^ + 1)(A^ + 2)/6 degrees of freedom. The 
total number of simulated maps must exceed N^ap > where jS^, is the 
variance of the skewness of a distribution. To estimate a skewness with vari- 
ance smaller than 0.1, the number of simulated maps has to be larger than 
6 X 10^ 



3 Topological Measures of a Random Field 

In this section, we will discuss the topological measure of a random field and 
the application to the test of deviation from Gaussianity. This approach was 
studied in |^||3^. It is found that among all these quantities, the Euler- 
Poincare (EP) characteristic is the most effective topological measure with 
regard to testing Gaussianity. Adler also derived the mean for a special 
non-Gaussian field: the field. Subsequently, Coles applied the result 
to a number of non-Gaussian fields which are derived from Gaussian. Gott 
et al. applied EP characteristic to simulated cosmic string maps. Both 



confirmed that the EP characteristic is effective in testing Gaussianity. We 
will briefly review the existing results obtained by previous investigators, 



then we will move on to study the statistics of the topological quantities. 
Special attention is paid to the real experimental situation where the noise 
level is high. We also derive some new results on topological measures in the 
presence of Gaussian noise. We will show that in the case where the signal 
to noise ratio is around 1:1, the EP characteristic, unfortunately, fails to be 
effective in discriminating between Gaussian and non-Gaussian distributions. 



3.1 Mean 

The central concept on which all topological measures are based is the excer- 
sion set ||2^ which is the set of points where the field F{r) exceeds a global 



value u. If we take a map of a certain area of the GBR sky, the excersion set 
of the map above level u, denoted by Su, will in general consist of a number 
of disjoint regions, each having a boundary which is the contour of satisfying 
F{r) = u. Since the smoothness and the differentiability of the contours are 
guaranteed by the beam smoothing^, the Euler-Poincare (EP) characteristic 
of Su, r„ is given by: 

Tu = ^ f kds, (22) 



2tt 

where k is the geodesic curvature of the contours and the intergal is taken 
over all contours in the map. 

For Gaussian random fields, the mean of the EP characteristic is exactly 
calculatable: 

< >= (^(^)''^^^P(-'^V2), (23) 

where Cq =< T'^{f) >, C2 =< V T{'r)'VT{f) > is the variance of the 
temperature anisotropy gradient^, and u = u/a, a is the rms temperature 
anisotropy. 

Among all possible non-Gaussian PDFs, we specifically choose the Xn 
distribution because of its relevance to the 0{N) cr- model. For our stan- 
dardized PDF given in Eq. (|TUp, the mean number of the EP characteristic 
is given by: 

<^^>- i^fl""^*'"""''^^ + 1) exp(-( Vi. + n)/2). 

(24) 



3.2 Statistics 

The Euler-Poincare characteristic is a discrete point process defined over the 
underlying random field. For any point process N{f) with mean < >= N, 

^For noise, without smoothing, the contour will be a fractal and it will not be dif- 
ferentiable. We illustrate this point in Fig. 2. For a detailed discussion, see citeadler, 
coles. 

•^C2/Co is equal to V2^'^/9l, where 7 and 6'* is the spectral parameters defined by p2| . 



the variance of the process is given by 



m 



< N'^ >=\N\+N\l + C), (25) 

where 

^=^1 dn,dn2ari2) (26) 

is the sample-area averaged two point correlation of the process A^. The first 
term in Eq. (^) is due to the discreteness of the process. The variance of 
the process is given by: 

=< > -N^ = N + N'^l (27) 

Thus, to find the variance of the EP characteristics, the central problem is 
to find the two point correlation. We approximate this correlation function 
^ by the peak-peak correlation x of the same underlying random field. Since 
r(z/) is defined through the z/-peaks of the density field, the approximation 
should provide a reasonable fit to the true F correlation. We expect that this 
approximation will give poor results when r(z/) ~ 0. 

The peak-peak correlation function is well studied in the literature |^ 
Q and is most easily understood as the following: Let Pi(f) be the proba- 



bility of finding one i/— peak at f and P2 be the probability of finding a u— 
peak at fi and f2. Then, the two point correlation for the peaks will be: 

l + C2{h,h)=P2/Pl (28) 
For a Gaussian random field, we have 



Pi = -^=^ I dxexpi- — ), (29) 



1 1*00 

7^ / dxexp{- — 
and 



27rdetM 



/ / dxidx2ex'p( — xM^^x), (30) 

J U(TQ J U(TQ 2 

where x = (xi,X2) and M(ri,r2) is the correlation matrix of the random 
field at ri,r2, which is 



M 



1 m 
m 1 



• <yo (31) 



and V'(^) is the normalized two point function. Thus, in the limit where 
ip <^1, the peak-peak correlation function is given by: 

ipeak = z/V(^)- (32) 

The beam-smoothed two point temperature correlation is well approximated 
by : 

C(0)=C(O)exp(-eV^c), (33) 



where the coherent angle ~ is a function of the beam width. The 
averaged correlation in an area A = 7i9\ is 



^ = Aj^/- (34) 

Thus, the variance of the Euler-Poincare characteristic is given by: 

a' = \N\ + u'i^rN',N = nin9\). (35) 

The mean and la uncertainty of the EP characteristic r(z/) for Gaussian 
random field is plotted in Fig. 3. Because of the approximations we used to 
derive the statistic, we expect the uncertainties of r(z/) given by Eq. (^) 
are not exact when ~ 0. 



3.3 The EP Characteristic of Noisy Data 

Our major concern over the applicability of the EP characteristic as a reliable 
statistic to discriminate between a Gaussian and a non-Gaussian random 
field is the noise term, which appears in all current GBR experiments and is 
comparable in amplitude to the signal. In this section, we will show first how 
the beam used in the experiments regulates the noise. Without smoothing, 
noise will be an obvious hazard to topological measures because the rms of 
the derivatives of noise is not well defined. We will then proceed to study 
the change of EP characteristic due to the beam-smoothed noise term. 
The spectral parameter of the noise is given by: 

= AY^{21 + 1) exp(-/(/ + l)^,^) ~ A{-)\ (36) 

1=2 '^s 

al = AY. l{l + 1)(2/ + 1) exp(-/(/ + ~ 1^(1)^ (37) 

1=2 ^ 

4 = ^E(^ - + m + 2)(2/ + 1) exp(-/(/ + l)al) ~ \A{-f. (38) 

1=2 

2 

Thus, 7 = = 0.75 and 6^ = = v^cr^. If we compare this value 

for noise with that of a Harrison- Zeldovich primordial perturbation, where 
7 ~ 0.5 and 9^: ~ 1.8a, we can conclude that the noise term regulated fairly 
well by the beam. 

However, even with a beam-smoothing, the EP characteristic of the ob- 
served random temperature anisotropy pattern is still changed due the pres- 
ence of noise. We first consider a situation where the sky is dominated by 
a quadrupole only, as we showed Fig. (4a). The EP characteristic is very 
simple for this map: F = 2 for u = —1, 1, 2. However, when one puts Gaus- 
sian noise of the same variance to the signal, the map is dominated by the 



feature of noise, which are shown in Fig. (4b). Even though we don't expect 
the temperature anisotropics to be a pure quadrupole, this example makes 
us cautious when using the EP characteristic. In the reahstic cases where 
the signal itself is also a random field which is non-Gaussian, the problem is 
harder. A full analysis of the EP characteristic for the sum of two general 
random fields is not tractable. Only one special case where both fields are 



Gaussian is solved |2^. In our problem, where one field is non-Gaussian, 
there is no ready-to-use result to apply. Thus, we solve this problem by the 
following strategy: the xL'^ = ••• is a class of distributions, ranging from 
very non-Gaussian (small n) to slightly deviating from Gaussian (large n), 
and the EP characteristic of this class of distribution is known. Thus, we 
first find the modified non-Gaussian PDF due to the presence of noise, then 
we find the best fit distribution to this modified PDF. The EP charac- 
teristic of the noisy temperature map is thus approximated to be the EP 
characteristic of the best fit distribution. The validity of this approach 
lies partly in the fact that if the difference between two PDFs goes to zero, 
the difference between the corresponding EP characteristic also goes to zero. 

We assume the noise is Gaussian and uncorrelated with the GBR anisotropy 
signals. The PDF for the noise is: 

Piv) = ^^expi-r^y2al) (39) 



where cr^ is the variance of noise. In the following, we consider only the case 
where the signal to noise ratio is 1:1, thus anoise = cr signal = l/"\/2- It is 
straightforward to generalize to the arbitrary noise case, and we will show it 
here. 

The PDF for the noisy map is given by: 

P{x) = J Psignai(a;')Pnoise(a; - x')dx' (40) 



It is convenient to use the cumulant function K(u), which is the logarithm 
of the characteristic function 0(u), the Fourier transformation of the PDF: 

(f){u) = J P(a;)e™^da;. (41) 

The cumulant function for P{x) is simply the sum of the cumulant function of 
the signal and the noise. For Gaussian noise, the cumulant function is simple: 
Knoise{u) = -^^^ signal, the PDF is the modified distribution. 

After some algebra, the cumulant is found to be: 

Ksignaliu) = - ln(l - ^) + (42) 



Thus, the cumulant function for the noisy signal is given by: 



Z XL 

K{U) = Kr,oise{u) + K^ignaliu) ^ -u" /2 + ——{ — )\ n > 1 (43) 



Thus, with noise, the cumulant is still that of a distribution in the limit 
of large n, but with degrees of freedom, where = ~ l.7n. In Fig. 
5, we plot the EP characteristic for n = 12 with and without noise. Even for 
the global monopole where n ~ 12, the EP characteristic fails to be effective. 



4 Three Point Temperature Correlation Func- 
tion 

There are various examples where non-Gaussian processes possess Gaussian 
PDF |T^. One classic example is the smoothed Poisson point process. The 
process is non-Gaussian when the smoothing scale is small, and tends to be 
Gaussian when the smoothing scale becomes large by the virtue of the central 



limit theorem [|T^. Thus, the test of Gaussianity should go beyond the mere 
one point PDF. The EP characteristic discussed in the previous section is 
one way to take into account the full properties of a random Gaussian field. 
But the noise present in the data prevents it from being effective. However, 



as we stressed before |]T5|, the three point temperature correlation function 
is a good measure of deviations from Gaussianity for the noisy data, as long 
as the noise is mutually independent and not correlated with signals. In this 
section, we will first introduce statistics to test these aspects of noise. Then, 
we will discuss theoretical predictions of the three point function in different 



models ||3J]. The reduced three-point functions for COBE 53 GHz signal and 
noise maps are obtained and show no deviation from Gaussian. We conclude 
that at the COBE scale, the temperature anisotropics are probably Gaussian. 



4.1 Properties of Noise 

A good understanding of noise in GBR experiments is crucial in testing the 
Gaussianity of the primordial density perturbation through the existing data. 
One has to make specific assumptions about the instrumental noise in order to 
test the Gaussianity: the instrumental noise must be mutually independent 
among pixels and Gaussian. The assumption has to be tested thoroughly 
before any attempt to decipher the non-Gaussian GBR signal from the data. 
The statistical tools of testing Gaussianity we introduced throughout this 
paper should also apply to the noise and we wouldn't repeat them here. In 
this section, we will address the following aspects related to the noise in the 
experimental data: (1) is the noise mutually independent? (2) does noise 
correlate with the GBR signals? 

In testing the mutual independence of the noise, we use the following 



results from statistics [25|: Let yi, {i = 1, A^) be a measurement of a zero 



mean random process. The mean and two point correlation of the process 



can be estimated as the following: 



1 ^ 



N 



Ey^^ (44) 



1=1 

I N-\s\ 

C2{s)^- ^ (|/,-y)(|/i+H-^),s = 0,±l;...,±(iV-l), (45) 



1=1 



p(s) = C2(s)/C2(0). (46) 



Then, if yt are mutually independent, then p{s) is asymptotically Gaussian. 
In particular, p(l) is asymptotically Gaussian with zero mean and variance 
j^. With a suitable redefinition of i/i, this result can be applied to answer 

questions (1) and (2). 

(1) Testing for the Independence of the Noise 

Let rii,i = 1, ...,N, be the noise in ?'th pixel. The mutual independence 
of noise can be tested through a second covariance of analysis of the square 
of the noise. Let i/i — rjf,i — 1, N. li rji are not correlated, then yj will 
not be either. Thus, we can define the following statistic: 



W^ = emm£, (47) 
'Ci(0, 1, 1) 



where 

1 AT-l 

CM i> 1) = ^ E ivi - yfivi^^ - yf- (48) 
i=i 

Under the hypothesis that the noise are independent, Wn is distributed as a 
standard Gaussian (zero mean and unity variance). We can define Wg, the 
same statistic as Wn by using the signal in each pixel, to test the independence 
of the signal. 

(2) Testing the Correlation among Signal and Noise 

Noise should not only be uncorrelated, but should also be independent of 
the signal. We assume that the noise is additive, i.e., 5obs = 5cbr + V- Once 
the noise is found to be independent, the variable 

yi-&'^-^rjl{5,-rm). (49) 

should also be mutually independent. Following the previous section, the 
following set of statistics 



pmcm^ (50) 

'C4(0, 1, 1) 



with 

^4(0, 1, 1) = ^ E iyi - y)'iyi+i - yf^ (^i) 



^ N-l 



is a standard Gaussian. 

(3) Analysis of COBE GHz DMR Map 

For COBE DMR, the (A+B) /2 is the signal+noise map and the (A-B) /2 
is the noise map. The properties of signal and noise are tested through the 
statistics defined aboved and shown in Table 2. 

The analysis shows that at the 1.6% confidence level, the noise is uncorre- 
lated. At the 19% confidence level, the noise is not correlated with the signal. 
Thus, we can conclude that the noise in COBE 53 GHz map is not correlated 
with the signal and marginally uncorrelated. The Wg in the last row is the 
statistic to test the independence of the signal. The statistic shows that the 
hypothesis that the data is not correlated failed badly. This is expected if 
the signal is primordial CBR fiuctuations. 

4.2 Three Point Function 

It is convenient to use the normalized two point and three point function 
7], where 

V^(|ri-f2|) = C2(fi,f2)/C2(0), 7/;(0) = l (52) 

V{fuh,h) =< ^(ri)^(r2)^(r3) > /CM'^ v{0) = f^s, (53) 

where /is is the skewness. The theoretical predictions for three point function 
are mainly the following: 

(1) Inflation 

Various non-standard infiation models ^ will generate a non-zero three 
point correlation function. The generic form of the three point function in 
most infiationary models is: 

Vih, r2, rg) = ^(^Adri - hD^Ufi - r^]) + 

i^i\ri - fiDipilh - ^sl) + i^ilh - fsDi^ilh - hi)), (54) 

where A is a dimensionless constant. As we will show in the Appendix, for one 
field slow-roll infiation models, A ~ 10~^. In non-standard infiation models 
[^], (7 can be much larger (up to order unity). 

(2) fields. 

Consider the field Y = J2i=i -^i- -^^ discussed before, the field 
describes the global topological defects in the large N limit. By extrapolating 
to low N, one can also gain insights into the possible temperature anisotropy 
patterns generated by domain walls, strings, monopoles or textures. The two 
point correlation is given by: 



^2(^,^2) = 2n0^(fi,f2), 



(55) 



where is the common covariance function of the Xj. We choose so that 
the two point function of Y matchs the observation, 

0'(ri,f2) = C2(ri,f2)/2n. (56) 

The three point function is found to be: 

'n{fi,r2,h) = \/8/n{ilj{\fi - hD-ipilh - hD'^^ilh - hDf^'^ (57) 
(3) Late-time Phase Transition 



In this model |3^, due to the conformal invariance of the system at the 
critical point [36], the three point function result from this class of cosmo- 
logical phase transitions has the following simple form: 

V{ri, rs) = a{tp{\ri - hD^Hh - r3\)^{\f2 - fsDf, (58) 

where a is a dimensionless constant of order unity. 

The full structure of t] is complicated. We consider the reduced three 
point function p7| where 



Vre,uce,{h,r2) =< ( (n) ^ (f^) > /^^(O)^/^. (59) 

The theoretical predictions in various models are given by: 

^redw(n,^2) = y(l + 2V^(|n - ^2!)) for inflationary models (60) 

Vreduced = ^/2/NiJ''i\h-r2\) for 0(N)or - modcls (61) 
'nreduced = aip^ilri-hl) for LTPT models (62) 

The three point functions from the previous three different categories are 
plotted in Fig. 6, the skewness is chosen to be the same for all cases. 

4.3 Analysis of COBE data 

The 53 GHz COBE DMR data is analyzed ultilizing the statistics we dis- 
cussed above. This frequency (53 GHz) is chosen because it has the best 
data quality [§. Before subtracting the dipole and beginning further analy- 
sis, the signal is weighted by the estimated pixel uncertainty. The dipole is 
subtracted by using the most recent COBE result |^. For 53 GHz (A+B)/2 
map, the subtracted dipole signal is: 

T{l,h){mK) = -0.198-cos(/)cos(6)-2.075-sin(/)cos(6)+2.333-sin(6), (63) 

where I and b are the galactic longitude and latitude. In our analysis, we 
consider only the 2019 pixels whose galactic latitude is 20 degrees or above. 
The COBE 53 GHz (A+B)/2 and (A-B)/2 are analyzed and the reduced 
three-point functions are shown in Fig. 7a, 7b. The result is consistent 
with Gaussian fluctuations and there is no deviation above the statistical 
uncertainty. We conclude that the current analysis shows that the statistics 
of CBR at the COBE scale are probably Gaussian. 



5 Statistical Tests on Intermediate Angular 
Scales {6 - 1°) 



The Gaussianity question is hard to resolve on the COBE scale {6 ~ 7°) even 
if there is no noise in the experimental data. One should expect the GBR on 
the large angular scale to be Gaussian simply by virtue of the central limit 



theorem ||T9[. Furthermore, there is intrinsic uncertainty in the statistical 



quantities measured in our local universe due to the cosmic variance 
which makes it harder to discriminate between Gaussian and non-Gaussian 
fluctuations on the COBE scale alone. On intermediate angular scales {6 ~ 
1°), the cosmic variance is small and the chance to detect a non-Gaussian 
signal is higher. As the data on these scales is accumulating, it is timely 
to consider seriously testing Gaussianity on these scales, since data in some 
experiments show clear non-Gaussian features. Although the current 



experiments are inclusive, due to the possible foreground contamination, it 
gives us hope that the Gaussianity problem will be resolved experimentally 
in the near future. The EP characteristic and the three point correlation 
will apply equally on both large and small angular scales if the sky coverage 
of the experiment is substantial. However, as the current state-of-the-art 
intermediate scale GBR experiments cover a tiny fraction of the sky, more 



sophisticated statistical tools [|IJ] are required to carry out the tests. In 
this section, we will introduce and discuss the bispectral analysis and the 
Hotelling's statistic. We will show that the statistic is a powerful 
statistical quantity to use on these scales and we also estimate the minimun 
data sample size to carry out the Gaussian test through statistic. 

In most current intermediate scale experiments, the data is sampled either 
in thin long strips or an annulus around an axis. In both cases, the data is 
one dimensional. The three-point function and bispectral analysis of 1-d 
data are well studied by statisticians and much of the specific techniques and 
mathematical details of this section are contained in monographs by p2[ 
which interested readers should consult. One can arrange the dataset as 
a time series, Xt, where in the present case, "time" t refer to successive 
positions in the sky. The data is usually edited so that the mean is removed. 
In this case, the three-point function is simply 

1 ^ 

Cti,t2 =< XtXt+tiXt+t2 >t= -TJ^^iXi+tiXi+ti, (64) 

i=l 

where the expression on the left-hand side is an estimate of three point func- 
tion from a data sample and is the size of the sample. Since the temper- 
ature anisotropies are always real-valued and assumed to be homogeneous 
and isotropic, the three point function has the following symmetry: 

6l,i2 = 62,*! = ^-tl,t2-ti = iti-t2-t2 (65) 



The bispectral density f{uJi, UJ2) is the Fourier transformation of ^t^.tj, where 



/(^l,^2) = I dhdt2e-''''^'-''''^'it,M. -TX < < 71. (66) 

In general, f{uJi, UJ2) is a complex function, following the symmetry of ^, one 
obtains the following symmetry relation for /: 

f{uji, UJ2) = fiuJi, Ui - UJ2) = fi-uJi - UJ2, UJ2) = f*i-uji, -UJ2) (67) 

Because of the symmetry, one just has to estimate the bispectral density in 
a small portion of {uji,u!2) parameter space. 

The unique feature of the Gaussian process is that the bispectral density 
vanishes for all u, i.e., /(ci;i,u;2) = for all ui,U2- To test this hypothesis, 
one can use the Hotelling's statistics p5|, which is constructed from the 



bispectral estimates defined on a "coarse-grained" frequency grid, {uji,ujj), 
where 

LJi = —,ujj = —,i = 1,2,..., L; j = i + l,...,j{i), (68) 

where L = [^], 7(0 = ~ [f] ~ 1- The parameter K is chosen to be much 
smaller than the sample size so that the frequency grid is "coarse-grained". 
Let rjij = f{(jji,ujj) and rearrange rjij into a vector rj = (xi,...,xp), where 
P = Z]iLi(7(0 ~ i) so that for each 1,1 < I < P, Xi = Vij- 

To estimate the bispectral density at each "coarse" grid point, one can 
construct a "fine" frequency grid around each {ui,ujj) point. Specifically, 

pDiT 

(^ip = uji + -J^^P = -'^^ -'^ + 1' •••''^ 
qDn 

uJi^ =uji + q = -r, -r + 1, r, (g ^ 0) (69) 

where the distance D is chosen so that the bispectral estimates at neighboring 
points on this fine grid are approximately uncorrelated. To insure that points 
in different "fine" grids don't overlap, it is required that D < j^^^^qifj- Since 
the total number of points in each "fine" grid is (4r + 1) and there are 
"coarse-grained" grids, thus the constraint on parameter r is: (4r + l)K^/3 < 
N. 

Let f{uip,Ui^) denote the estimated bispectral density function at the 
points (co'jp,ci;jJ. Due to the careful choice of grid point, one may regard 
the set of estimators {/(tUj^, cjj^)} as n = 4r + 1 uncorrelated and unbiased 
estimates of f{u!i,ujj). Forming the bispectral estimates f{ijji^,u}i^) into a n 
column vector, denoted by ^ = (^1,^2, •••,^n), then, at each "coarse-grained" 
grid point xu there will be an estimated bispectral density ■* from the "fine- 
grid". When the sample size N is large, = 1,2, ...,n), is distributed as 



complex normal with mean rj and covariance matrix S^. The maximum 
likelihood estimates of rj and are: 



-in A n 

^ = -Ed^ S = -, A,, = 5:((6,-r/)(e,,-r^r). (70) 
^ i=i ^ fc=i 

The Hotelling's statistics are defined as: 

= nr/A^^r/. (71) 
Under the Gaussian assumption, the mean vector r] = and the statistic 

F = ^^^T^ (72) 

is distributed as a central F-distribution with (2P, 2(n-P)) degrees of freedom. 
To test Gaussianity through degree scale experiments, the sample size 



should be large enough to avoid "sample variance" |^3| . The minimal sample 
size that could be used to carry out the Gaussian statistical tests can be 
estimated as follows: to carry out the bispectral analysis and Hotelling's 
statistics, one should use at least two "coarse" grids, and because the cor- 
relation angle of the beam-smoothed GBR anisotropies is around 6^,, where 



1.86s for Harrison- Zeldovich primordial spectrum |T2[. Thus, we need 



at least two "fine" grids on either side of the "coarse" grid. The minimum 
parameters to test Gaussianity is L > 2 and r > 2, which in turn gives 
K > 4. Thus, the minimum sample size is: > (4r -|- l)K'^/3 = 48. This 



is about twice the current largest sample at the half degree scale This 
result is encouraging: it suggests that we don't need a full sky coverage with 
half-degree resolution to carry out the Gaussianity test, and the Gaussian- 
ity question can be resolved by ground-based experiments or balloon fiights. 
With regard to sampling in the sky, the data should sample as sparsely as 
possible to avoid possible correlation between different data points, and it 
will be better to sample in thin long stripes or annuli to take advantage 
of the simplicity of one dimensional dataset. Currently, most degree scale 
experiments have fewer data points and are thus inadequate to perform bis- 
pectral analysis. However, we expect the situation will change dramatically 
soon, and the results we discuss here will be helpful for designing future 
experiments. 



6 Gaussianity of CBR on Small Scales 

In this section. We will show how to use the statistic by analyzing the 
recent 96 point RING data from OVRO . A special class of non-Gaussian 
signal, point sources, which are interesting in their own right, is separated 
out and discussed in detail in the next section. The data sample is weighted 



according to the error in each pixel. The sample size N = 96. The esti- 
mate of the three-point function is shown is Fig. 8a and the the real and 
imaginary part the bispectrum is shown in Fig. 8b, 8c. We have used an 



optimum window function (see ||2^ for details) to smooth the discrete data. 



For Gaussian distributed data, both the three-point function and bispec- 
trum should vanish. The results shown in Fig. 8a, 8b, 8c already suggest 
that the data might be non-Gaussianly distributed. To show how statisti- 
cally significant the deviation from Gaussian distribution, we can use the 
Hotelling's statistics. The parameters we choose for RING data are: 
K = 4:,r = 2,n = 9, d = 4, P = 2, L = 2. With this choice of parameters, 
the statistical distribution for F is shown in Fig. 9. The 95% confidence 
level (C.L.) upper limit of F is Fc = 3.15 if the data are Gaussian dis- 
tributed. Hotelling's statistic estimated from RING data is = 5.76, or 
F = 7T^/2 = 20.1 which is much larger than the 95% C.L. upper limits for 
Gaussian distributions. Thus, we conclude that the data are probably not 
consistent with Gaussian statistics. 

Even though the data failed the statistical tests we proposed, one can- 
not conclude that the non-Gaussianity is due to non-trivial interactions in 
the inflationary cosmology or the topological defects produced in the early 
universe. The non-linear gravitational evolution will produce non-Gaussian 
signals which have to be carefully studied and subtracted to gain some knowl- 
edge about the Gaussianity of primordial perturbations. Part of the answer 
to this important issue is contained in the next section and we will not discuss 
here. 



7 Point Sources 

On intermediate angular scales (~ 1°), the current datasets are too small to 
carry out the three-point correlation and the angular bispectrum analyses 
we developed in the previous sections. To test Gaussianity based on the 
small dataset available (usually about 10 data points), the statistics have to 



be very custom-designed to be useful |T1]. However, a clear non-Gaussian 



signature will be the point-like GBR anisotropies. In fact, two candidates of 
such sources are detected in the MSAM experiment |^ . One source, located 



in a dust-free region, has a flux of 3.7±0.9Jy at 5.6cm~^. Another candidate 
has a flux of 2.9 ± 0.7Jy, is located 4.3° away from the flrst sources. Both 
sources are compact and have angular size less than the beam width a = 
0A259fwhm = 12 . Assuming that the angular size of the sources are half 
the beam width, one can flnd the the flux intensities are (3.8±0.8) x 10^ JySr~^ 
and (3 ± 0.6) x lO^JySr^^. As a comparison, the flux of a 6T temperature 
fluctuation will produce a flux 

exp(Ju^ikTo) hu ST 
explhu / kTo) — 1 kTq Tq 



where Tq = 2.73K is the CBR temperature and B,y is the CBR flux at fre- 
quency u. At 5.6 cm~^, the CBR flux is = 1.5 x lO^JySr"^. Thus, for a 
temperature anisotropy of AO^K, which is the theoretically estimated tem- 
perature anisotropy at the half degree angular scales for CDM with standard 



recombination the expected flux is 



h = 6.8 X lO^(-i^)JySr-i. (74) 

Thus, the first sources correspond to a 5a and the second sources correspond 
to a 4cr peak. 

In the Gaussian picture, for a z/— peak, the mean size 6 and the average 
distance between peaks d are given by : 




e ^ y/2—Jl- for z/>l, (75) 
71/ V 

d = 2/3(7rn^)-^/^ (76) 

where 9^: = 1. 29s and n^, is the number density of u peaks. For large u, riy is 
given by: 

^^ = ^%^«-P(-V2). (77) 

Thus, for 4(T peaks, the mean size will be 0.716'^,, which is marginally consis- 
tent with observation, but the the mean distance between rare peaks will be: 
~ 50^s ~ 20°, which is much larger than the angular separation between the 
sources. One may try to explain both point sources by Gaussian statistics, 
assuming that they are just 3 a peaks and fit the low limit of the observed 
flux. Then the average distance between peaks is 119 s-, which is roughly the 
same as the observed value. However, in this case, the averaged angular sizes 
of the peaks are 1.2^^, which is larger than the beam width. We conclude 
that if these sources are of CBR origin, they are not consistent with Gaussian 
statistics. 



Various topological defects, notably soft domain wall bubbles [^, the 



global monopoles or textures |]2^, are capable of producing spot-like 
CBR anisotropics of any size by appropriately choosing model parameters. 
However, before one rehes on topological defects as an answer, one has to 
filter out the foreground contaminations carefully. Several types of radiation 
may contribute the point-like sources observed in the experiments. One of 
them is the Sunyaev-Zeldovich (SZ) effect from rich clusters. The scattering 
of microwave photons by hot electrons in the intracluster gas will make a 
cluster a powerful source of submillimeter radiation. The typical angular 
size of the core of the hot gas is of order arcminutes, and the flux density is 
given by fg^: 

Fu = - 4) ■ ^ ■ S.(To), (78) 



where x = and y = J -^^aTUedl. For a typical rich cluster, Te ~ 
lO^i^', rig ~ 10~'^cm~^ and / ~ IMpc, the estimated y parameter is around 
10~^ — 10~^, The flux density is around 3 x lO^JySr"^ which is exactly what 
the MSAM experiment observed. Thus, it is very likely that the observed 
sources are due to SZ effects of unresolved rich clusters in the field. 

The multifrequency channel methods is widely used to separate the local 
contamination from true CBR signals. The method may not be effective to 
single out the SZ effect. As we showed in Fig. 10, the spectral index of the 
SZ effect is very close to that of primordial temperature fluctuations at low 
frequencies (z/ < 100 GHz). When the frequency gets higher (z/ > 200 GHz), 
there is a small deviation, but at this frequency range, the dust emission 
will dominate. A possible way to discriminate the SZ effect is through po- 
larization of the radiation from candidate point sources. Due to the peculiar 
velocity of the rich clusters, the radiation will be polarized in SZ effect 



However, the point-like CBR anisotropies from topological defects will not. 
But since it is currently hard to measure the polarization of the radiation 
down to required accuracy, we won't discuss this approach in detail here. 

Apart from the SZ effect, primeval dust or a population of IRAS- 
like galaxies at high redshift |^ may have substantial contributions at the 



submillimeter range. For dust grains, if one assumes the emissivity of the 
dust is ~ I/", a 1.5, the flux spectrum of dust emission is given by: 

/. - {^r^--^:^- (79) 

The peak of the distribution is located at 

Vpk = 4:.5kTd/h = 3750GHz(rrf/40K). (80) 

Thus, in order for the peak of a hot dust spectrum (T^ = 407^^) to be red- 
shifted into the 300 GHz range which the MSAM experiments operates, the 
redshift of the epoch of formation of primeval galaxies should be around 
z 10. The angular size of the dust envelope is 

0=^-^^, D = 2H,\1-^), (81) 

V 1 + ^ 

where I is the proper size of the dust envelope, which is about 10 — 100 kpc. 
Thus the typical angular size of the possible point sources produced by the 
primeval galaxies is 

e = 0.03 — — — ?. (82) 

lOkpc 10 ^0.5^ ^ ^ 

The observed flux density is 



(83) 



Assuming that, in the rest frame of the sources, the intrinsic flux is peaked 
around 3000 GHz and the luminosity is L, then the flux density is 



.r. 3t n 1/ L 10 3000GHz, 
Both the angular size and the flux density are too small to account for the 



observed flux in the MSAM experiment Thus we can conclude that 

primeval dust or the distant infrared galaxies do not account for the point- 
like sources observed. 

As we showed in Fig. 11, which is plotted according to Eq. (^), it is clear 
that if there is a population of cold dust (T^ ~ 4i^), then the flux density 
will peak near where the experiment operates. A uniform background of 
such cold dust is impossible unless the density the low enough so that the 
optical light from distant quasars will not be absorbed. However, clumpy 
cold dust is helpful to explain the experimentally observed point-like sources. 
A possible scenario to explain the spot-sources based on cold dust is the 
following: there is a population of very quiet galaxies where most of the star 
formation activities are shut down, so that there will be very low intensity 
radiation in the far-infrared regime. Thus, this population of galaxies is 
not observed by the IRAS flux limited survey. However, as the hot dust 
cools down to around 4K, they become powerful submillimeter emitters. The 
spectrum of cold dust emission is shown Fig. 11. Multi-spectral analyses can 
shed light on the possible spectral parameter of the sources. The problem 
with the cold dust scenario is that as one look back in time, these sources 
used to be very powerful infrared sources because the radiated flux oc T^^". 
Due to the abundance one observes today, we can estimate the luminosity 
at high redshift, which is much brighter. However, the newest result from 
the COBE FIRAS has already put a stringent limit on the possible 
evolution of infrared galaxy luminosity function. Detailed modelling is in 
progress. But we are pessimistic about explaining the point-sources based 
on cold dust scenarios. 

In conclusion, the point-like GBR anisotropy is a clear non-Gaussian sig- 
nature. If the future studies confirm that the point sources candidates are 
truly GBR fluctuations, it will be an exciting new chapter. It will provide 
direct evidence that topological remnants left over from the early universe 
do exist. 



8 Conclusion 

In summary, the following points which are related to the test of the Gaussian 
nature of the primordial fluctuations are discussed in this paper: 
(1) We listed the skewness and kurtosis in various physically motivated mod- 
els, with and without noise. We also discussed the use of the multivariant 



skewness and kurtosis to quantify the deviation of a distribution from Gaus- 
sian. 

(2) We discussed in detail the Euler-Poincare characteristic of random fields. 
We showed that the Euler-Poincare characteristic will not be a good discrim- 
inator between Gaussian and non-Gaussian random field when the noise is 
comparable to the signal. 

(3) We stress the use of the three-point temperature correlation function to 
test Gaussianity. The predictions from various models are discussed and the 
COBE 53 Ghz data is analyzed. The analysis shows that the fluctuations 

are probably Gaussian on the COBE scale. 

(4) We discussed the detailed statistical tests on intermediate angular scales. 
The bispectral analysis and Hotelling's statistics are emphasized. We also 
discussed briefly the sampling technique and minimum sample size to test 
Gaussianity statistically on half degree scales. 

(5) We discussed testing Gaussianity on small angular scales (arcminute 
scales). The RING data from OVRO is analyzed and shown to be propably 
non-Gaussian. 

(6) We discussed looking for point-like sources as a way to test Gaussianity 
and hunt for topological defects in small scale GBR experiments. The SZ 

effects and the effects of primeval dusts are discussed. 

The current status of testing Gaussian nature of GBR anisotropics are 
summarized in Table 3. 
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Appendix 



A Gaussian Nature of Perturbations from In- 
flation 

In this appendix, we will discuss the Gaussianity of the primordial fluctu- 
ations produced in inflation. We pick one simple model, Linde's chaotic 
inflation model |ET| to analyze. The approach is rather heuristic but brings 



into focus the Gaussianity problem. For more rigorous treatment, see ^2 . 

The basics of the inflationary dynamics are the following: there exists an 
epoch where the universe is dominated by the vacuum energy of a scalar field 
(J). The Friedman equation which describes the evolution of the background 
metric is: 

H' = ^-^p, p = 1/20^ + ^^(0), (85) 
and the dynamical evolution of field is: 

^ + 3/f^+^ = 0, (86) 

where the inflaton potential for chaotic inflation is simply 

V{<f>) = X<j>\ (87) 

The model is easy to analyze in the slow-roll regions where ^ iJ0, which 
is sastified for > nipi/^/Tn. By using the gauge- invariant parameter ( 
introduced by where 



4 « ^m, (88) 

A0 is the zero-point quantum fluctuation of the field along its classical 
trajectory. The probability density function for A0 is a Gaussian with zero 
mean and variance given by (A0)^ = (^)^- Let us denote 2z^^ by then x 
is a random variable with standard normal distribution. 
Expressed in terms of x, ( is given by: 

Cfc = ■ X (89) 



In the slow approximation, the quantum fluctuation of is negligible. Thus, 
if the Hubble parameter is a constant during inflation (which is the case for 
original exponential inflation), one can conclude that the primordial den- 
sity fluatuation is Gaussian. However, since H is related to locally, the 
fluatuation in will give rise to a fluctuation in H. Taking this into account. 



where He is the classical value and nipi = is the Planck mass. The 
usefulness of the ( parameter lies in the fact that it is a constant throughout 
the inflationary, radiation and matter dominated epoch, and the Sachs- Wolfe 
contribution to the CBR temperature fluctuation is given by ||5^ : 



^I.i (91) 

T 5 ^ ^ 



Thus, the large angular temperature anisotropy is given by: 

6T 




T 50 ^ ' V 127r rripi ' ''^'^^ 

The statistics of CBR fluctuations are non-Gaussian because of the term. 
This expression for the deviation from Gaussianity is rather generic in in- 
flationary models and will lead to the functional form of the three point 
correlation we will give in section 4. The difference is that various models 
will produce different skewness. In the one-field chaotic inflation model we 
treat now, the skewness is given by: 




(93) 

127r rUpi 



During the slow roll [17 



^ = aJ2/3X'/^N!/\ (94) 


where ~ 60 is the number of e— folds the scale factor inflates during 
inflation. In order for the amplitude of the fluctuation to be of the same order 
of magnitude as the CORE observation [Q, the self-coupling constant A of 
the inflaton field is given by A ~ 10~^^. For such a weakly coupled field, the 
deviation from Gaussian is estimated to be: /is ~ 10^'' with ~ 0j ~ 5mpi, 
which is negligibly small. 

In conclusion, we have shown that the primordial fluctuation from slow 
roll inflation is very close to a Gaussian. The deviation from Gaussianity is 
of order 10~^ in the chaotic inflation model, and the result is general to other 
one-field inflation models. The root cause is the small coupling constant 
of the inflaton potential, which is required to give the right amplitude for 
primordial fluctuations. Thus, Gaussian-distributed CBR anisotropics are a 
nature result of lots of inflation models. 
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Table Captions 



Table 1: The skewness and kurtosis of the temperature anisotropies in various 
non-Gaussian models. 

Table 2: Testing the properites of signal and noise in COBE 53 GHz map. 
Wn is the statistic for testing mutual independence of noise. Under the null 
hypothesis that the noise is uncorrected, Wn will distribute as a standard 
normal. Wc is for testing correlation between signal and noise, and Wg is 
for testing the correlation between signals. Under the null hypothesis (signal 
and noise are not correlated; signals are not correlated), both Wc and Wg are 
distributed as standard normals. 

Table 3: A summary of current status on testing Gaussian nature of GBR 
temperature anisotropies on different angular scales. 



Table 1 



Non-Gaussian Sources 


No Noise 


Signal to Noise Ratio 1:1 


Skewness 


Kurtosis 


Skewness 


Kurtosis 


Soft Domain Wall^ 


1.4 


3.75 


0.5 


0.94 


Cosmic string^ 





1.5 





0.38 


Global monopole"^ 


0.82 


1.25 


0.29 


0.31 


Global Texture'^ 


0.71 


0.94 


0.25 


0.25 


0{N) (T-modeP 




15/4iV 


^l/4iV 


15/16iV 


SZ from Rich Cluster^ 


-0.66 


2.72 


-.28 


0.22 



^Extrapolated from 0{N) model 

^On Scales of several arcminutes 

■^COBE sensitive scale, no beam-smoothing 

''On Scales of several arcseconds to arcminutes 



Table 2 





estimate 


Confidence Level 




2.54 


1.58 X 10-2 




1.21 


0.19 




8.73 


1.1 X 10"^'' 



Table 3 



Angular scale 


Small 


Intermediate 


Large 




(~ arcminutes) 


(~ 1°) 


(> 2°) 






South Pole, 




Current 




MAX, MSAM, 


COBE 


Experiments 


OVRO 


SK93, PATHON, 
White Dish,... 


MIT balloon 






Inflation: 




Theoretical 




Gaussian; 




expectations 


Non- Gaussian 


Defects: 
non-Gaussian. 


Gaussian 


Analysis 




Sky coverage 




from 


Non-Gaussian 


is still 


Gaussian 


experiments 




too small 






Source 


Eventually 


Central limit theorem; 


Comments 


contaminations 


Decisive 


Cosmic variance 



Figure Captions 

Fig. 1: The probability distribution functions of the modified distribution 
with n degrees of freedom. The dotted hne is for n=4, the short dash hne is 
for n=8, the lone dash line is for n = 16 and the solid line is the standard 
normal distribution. The values for n are chosen because the corresponding 
distributions describe the domain wall, global string and texture, respectively. 
Fig. 2: The simulated unsmoothed Gaussian noise 5 in a 10° x 10° patch. 
The solid line is for the contour 5 = 0, the dotted line is for the contour 5 = a 
and the short-dashed line is for 5 = 2cr, where a is the standard deviation of 
the noise. 

Fig. 3: The statistic of EP characteristic for a Gaussian random field. The 
soild line is the mean EP characteristic and the dotted line is the la uncer- 
tainty estimated from Eq. (^). The unit of the vertical axis is arbitary. 
Fig. 4a: The contour plot of the sky signals S if they are dominated by a 
quadrupole. The solid line is for the contour S = 0, the dotted line is for the 
contour 6 = ±0" and the short-dashed line is for 6 = ±2ct, where a is the 
standard deviation of the observed sky signals. 

Fig. 4b: The contour plot of the sky signal 6 when the signal to noise ratio is 
1:1, where the signal is a quadrupole. The solid line is for the contour 6 = 0, 
the dotted line is for the contour S = ±a and the short-dashed line is for 
6 = ±2a, where cr is the standard deviation of the observed sky signals. 
Fig. 5: The genus curve for distributions with n degrees of freedom. The 
solid line is for the random Gaussian field, the dotted line is for n = 12 and 
the dashed line is for n = 21, which describes the noisy case. 
Fig. 6: The theoretical prediction of reduced three-point function in three 
different models: The dotted line is for inflation, the short dash line is for 
0{N) a model and the long dash line is for LTPT. The reduced two point 
function, which is the solid line, is modeled as exp(— 6*^/6'^), where 6c = 13.5°. 
Fig. 7a: The reduced three point-functions estimated from COBE DMR 53 
GHz (A+B)/2 map and (A-B)/2 map. The sohd line is for the (A+B)/2 
map and the dotted line is for the (A-B)/2 map. 

Fig. 7b: The estimate and error of the three point function for COBE 53 
GHz (A+B)/2 map. The result is consistent with prediction from a Gaussian 
field (which is zero). 

Fig. 8a: The reduced three-point function estimated from RING data. We 
place 96 RING data points on a one dimensional lattice and the horizontal 
axis is the number of lattice spacing between the points used for estimating 
the reduced three point function. 

Fig. 8b: The plot of the real part of the estimated bispectrum f{uJi,uj2) 
estimated from the RING data. In this plot, we choose Ui = 002 = w. The 
frequencies u is plotted in unit of 7r/20. 

Fig. 8c: The imaginary part of the estimated bispectrum f{ui,LJ2) from the 
RING data. In this plot, we choose uji = LJ2 = The frequency uj is plotted 
in units of 7r/20. 



Fig. 9: The statistical distribution oi F = ^^^p^T^. We choose P = 2 and 
n = 9 for the RING data. 

Fig. 10: The frequency dependence of the antenna temperature for CBR and 
for SZ effect. The soild hne is for CBR and the dotted hne is for SZ. The 

antenna temperature is normahzed so that it is unity at low frequencies. 
Fig. 11: The flux density of the 4K cold dust. The cmissivity of the dust is 
chosen to be ~ i^^'^. The unit of the vertical axis is arbitary. 



